#include "saxpy.h"

#include <arm_sve.h>

#include <stdexcept>

void saxpy(const float* x, float* y, float a, const size_t n) {
  for (size_t i = 0; i < n; i++) {
    y[i] = a * x[i] + y[i];
  }
}

void saxpy_sve(const float* x, float* y, float a, const size_t n) {
  // SVE implementation is not yet complete.
  throw std::runtime_error("SVE saxpy implementation is not yet complete.");
}